The magneto-optical Barnett effect and spin momentum transfer 
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The interaction of polarized light with a spin in the presence of dissipation is shown to be equiva- 
lent to a spin transfer process that can cause switching. In plasmas, the spin transfer is dominated by 
a spin-spin exchange term while at lower energy densities it is dominated by an optical Barnett-like 
effect. This latter effect is used in conjunction with optical phonons to predict femtosecond magne- 
tization reversal believed to be recently measured in GdCoFe thin films. Conventional approaches 
based on the Bloch and the Landau-Lifshitz equations do not reproduce this ultrafast switching. 

PACS numbers: 76.60.+q,78.20.Ls,75.40.Gb, 76.60.Es, 52.38.-r 



In 1908, Richardson proposed of what has come to be 
known as the Einstein-de Haas effect By conserva- 
tion of angular momentum, a body which is free to ro- 
tate around an axis will indeed start rotating upon mag- 
netization. A few years later, the inverse effect, that is 
magnetization by mechanical rotation, was proposed and 
observed by Barnett Q ■ Barnett showed that such rota- 
tion is equivalent to a magnetic field directed along the 
axis of rotation and effectively contributes to the inter- 
nal field in a magnet. The direction of magnetization in 
a body at rest was shown to differ from that in a rotating 
body. A similar effect was also recently observed in para- 
magnets subjected to rotating magnetic fields instead of 
mechanical rotation Rotating black holes also emit 
particles (Hawking radiation) with net polarization due 
to the angular momentum of the black hole and can also 
be understood as a manifestation of a gravitational Bar- 
nett effect [4(. 

In this paper, we show that there is an optical version 
of the Barnett effect that can lead to magnetization re- 
versal in the femtosecond regime. By optical we mean 
that the laser imparts spin, rather than orbital momen- 
tum, to the magnetization. The optical Barnett effect is 
shown to act in similar fashion to a direct spin momen- 
tum transfer between the laser and the spin polarization 
of a plasma at high power. 

We believe that the optical Barnett effect has already 
been observed in a recent experiment [B| . The experiment 
demonstrated that it is possible to switch the magnetiza- 
tion of a GdCoFe film with a circularly polarized, 40 fs, 
800 nm laser pulse at relatively low power, of the order of 
10 11 W/cm 2 . Using a single spin picture, we show rigor- 
ously that the optical Barnett effect in combination with 
coupling to an appropriate bath, can provide a consistent 
explanation of the observed femtosecond switching result. 
Since in this case a real understanding of the short time 
response of the magnetization to the laser is crucial, we 
refrain from using the phenomenological modified Bloch 
equations (MB) [3] or the Landau-Lifshitz (LL) equa- 
tion but instead use a self consistent formulation for the 
spin and its environment. 

First we provide a motivation for the Hamiltonian that 



will be adopted for the interaction of the laser with the 
electronic degrees of freedom. Starting from the Dirac 
equation, Foldy and Wouthuysen Q used a series of 
canonical transformations to derive the Pauli equation 
and other terms of higher order in inverse mass. To or- 
der 1/toq, the Hamiltonian is given by 
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This Hamiltonian can also be derived using the gauge- 
invariant proper-time method js}]. The terms in this ex- 
pansion are well known. The third and fourth are the 
Zeeman and spin-orbit coupling energies. The fifth, less 
well known, term is a Heisenberg exchange-like interac- 
tion between the spin of the laser pulse and that of the 
electron. This term becomes important at laser powers 
of 10 15 W/cm 2 and above. 

A polarized laser pulse carries angular momentum 
but more importantly it also carries spin angular momen- 
tum. The total angular momentum J of the laser beam 
occupying volume V is proportional to the momentum 
E x B of the wave and given by (cgs) 



4ttc J, 



d 3 r rx (E x B) 
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This vector can be decomposed into two parts according 
to J = L + S, with L = 1/4ttc J v d 3 rE • (r x V) A being 
the orbital contribution, while the spin contribution is 



S = 



4-7TC 



d?r E x A. 



(3) 



For circularly polarized light, with polarization vectors 
e± = l/v2(x ± iy), and a vector potential A(t, r) = 
J d 3 k/(2n) {ek±a± (k) exp [i (k- r — u>t)] + c.c}, the 
spin angular momentum can be written as 

S = J- /'^k[|a + (k)| 2 -| a _(k)| 2 l . (4) 
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For a laser with B/ = (Bi cos , Bi smujit, 0), a_(k) = 
and a + (k) = correspond to positive and negative fre- 
quencies u>i , respectively. This form clearly shows the 
chirality of the laser wave and most importantly that S 
is parallel to the k- vector which is perpendicular to the 
B field. For completeness, we also point out that for 
beams with a finite cross section, the laser will have an 
additional component of B in the direction of propaga- 
tion of the wave which is proportional to the gradient of 
the transverse components [9j. In the following we will 
ignore this last contribution since it is important only at 
the edges of the beam. Moreover the orbital degrees of 
freedom will be considered as part of the environment. 

In order to give a clear discussion of the optical Bar- 
nett effect, we introduce an effective Hamiltonian, 7i cff , 
for the spin degrees of freedom derived from Eq.[T] TL e ^ 
includes energy exchange between three different sub- 
systems (cf. Fig.[TJ inset) and is {% = 1) 

H eff =H S +H„+Hs 
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The spin Hamiltonian, H s — — 7B; • a — ^Aa 2 zl in- 
cludes the interaction with the laser field B/ and an axial 
anisotropy term. The cr's are Pauli spin matrices for spin 
1/2. The anisotropy term is taken in the mean field ap- 
proximation to avoid having it trivial for spin 1/2 fields. 
It is only relevant after the the laser source is off. The 

gyromagnetic ratio 7 is assumed positive. The Hamilto- 

2 2 

nian, H q = ^- + ^-q 2 , represents a single optical phonon 
mode with energy w II Oil ■ The Hamiltonian TLq is that of 
the macroscopic bath [111 ]. The interaction Hamiltonian 
T~(.sqi — — ^S-a — X(E)a q — q Q, includes an exchange 
term and a linear coupling to the optical phonon mode. 
The effect of the electric field on the mode q is implicit in 
this interaction. The coupling constant in the exchange 
term is 



J = 



(6) 



This constant is very small but since |S| is proportional 
to the power, JS gives rise to fields > 1 Tesla for pow- 
ers > 10 15 W/cm 2 . While this exchange is negligible for 
the experimental conditions in Rcf. [5(, it becomes im- 



portant in plasmas 12]. The interaction between the 



spin and the bath is mediated by the spin-orbit coupling. 
The material-dependent parameter X(E) is proportional 
to the power of the laser which is responsible for the ex- 
citation^ the optical modes [13] . As first pointed out 



in Ref . [13J > a strong spin-orbit coupling can provide a 



fast relaxation channel for the spin of the electrons in 
the femtosecond regime. 

The coupled spin-laser-bath system is better studied 
in a frame (xi,X2,X3 = z) rotating around the z— axis 
with frequency lui. In this frame, the Larmor torque is 
time-independent and the spin Hamiltonian becomes 
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FIG. 1: Dynamics of M z , the component of the magneti- 
zation M along the Barnett field H_g, induced by circularly 
polarized light with frequency u>i. For curves cl, c2, c3 and 
c4, lui = W 15 Hz . For c3-, Ui = -W 15 Hz. The laser field 
B t is 10" 3 i/ s for c3 and c3-, 2 x 10" 3 i? s for c4, 1Q~ 5 H B 
for cl, and 5 x 10~ 4 //g for c2. M z switches from — 1(+1) 
to +1( — 1) for positive (negative) uoi as shown by the curves 
c3, c3— . The curve c3n shows no switching occurs for uii and 
M z (0) both positive. The different energy transfer channels 
involved in the excitation of the magnetization M by the laser 
uii and its fast relaxation by the macroscopic bath Q via an 
optical mode q with energy u>o ~ are shown in the inset. 
Presented results are for ujo = 0.8wi, T = 0.2, and A = 10 4 Oe. 



where Hg = uji/jz will be called the Barnett field and 
B = A(a z ) + JS Z . The equations of motion for (M), 
the average of the spin operator, are 

(AIi) = -7OB0 + H B ){M 2 ) + 2\{{M 2 h) - (M362)) (8) 

(M 2 ) - 7 (B + H B )(M X ) - 7 Bi<M 3 ) + 2A((M 3 M - (Ahb 3 )) 

(M 3 ) = 7 B,<Ma) - 2A «M a 6i) - (A/162)) , 

where b^s are the harmonic mode variables in the rotat- 
ing frame. The equations of motion for q are 



(^+u 2 )(q i ) = (Qi)+\(E)(a i ) 



(9) 
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M 1 -\Z I M l b l (t)-^-(B + H B )M 3 (7) 



The equations for the average qi are solved by treating 
the spin source as a perturbation. The bath Q is assumed 
Ohmic and is the source of dissipation in the mode q 
which we take to be T = 0.2 using Shcn and Bloembergen 
notation [loj . In the adiabatic limit and in the absence 
of anisotropy, the system spin plus bath, H s + 'Hq, has 
been treated earlier, see e.g. Ref. [l4(. Here we study the 
more difficult and experimentally relevant case of non- 
adiabatic magnetization reversal. 
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In the rotating frame, the instantaneous torque is 
modified by the Barnett field. For visible light with 
uji w 10 15 Hz, the optical Barnett field is as large as 
« 10 7 Oe. This is larger than most exchange fields and 
hence a single spin picture should be adequate even for 
the treatment of the interaction of a laser with a ferro- 
magnet. Despite its tremendous magnitude, the Barnett 
field does not induce femtosecond magnetization reversal 
unless three key requirements are met. First, the Barnett 
field has to be much larger than the laser field B;. Sec- 
ondly, the magnetization has to be coupled to at least one 
optical mode q of energy ojq w jHb ■ Third, the damping 
of this mode due to its interaction with the macroscopic 
bath Q has to enable efficient energy transfer from the 
spins to the mode q. All three requirements can be met 
by various combinations of the parameters, Hb, Bi, ujq 
and A(-E), but the range of each individual parameter 
depends on the particular choice of the others. 

Instead of going to the rotating frame, it is also possible 
to go to a frame which is rotating around the axis of the 
effective field H e ff = B; + (Hb + B ) z [HI]. In this case 
the transformation is explicitly dependent on the Barnett 
field and is given by 



U = U y (9) U z (w.) Uy 1 (9) U z (ui) , 



(10) 



with U a (e) = exp (— ieta a /2), u> s — jH e fi and 6 = 
tan -1 ( Bo ^} Hb J ■ It can be shown that this transfor- 
mation IA is a spin gauge transformation of the full 
U(l) x SU(2) symmetry of the Hamiltonian in Eq. [1] 

m 

To solve equations [8] and [9J we need to calculate the 
average of the product of two operators (&i(t)qj(t')) . For 
this we need the density matrix of the whole system or we 
may use instead the more powerful functional formalism 
(for a detailed discussion of this method and its appli- 
cation to sd exchange in metals see [l7| and references 
therein). The generating functional is 

Z[3x,3 2 ] = I^J DrjDpDqexp (^-i J dt (H s {f),t\l\) 
+H q (p, q) + Hsqi - Ji • q - J 2 • o))) Q 

where the spin variables are written in terms of Grass- 
mann variables, a = — |ry x 77 which allows the use of 
Wick's theorem in the path integral expansion and 3i 
are two virtual external sources. The average values are 
found in the usual way 17, 18| 
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Before switching on the laser, the phonon is assumed to 
be in equilibrium and is not coupled to the spins. After 
that, the spin is driven out of equilibrium by the laser 
and q is treated as a perturbation which responds lin- 
early to any changes in the spin. As in the sd exchange 



problem [l7|, the equations of motion are non-local in 
time and integrating out the optical mode q gives rise to 
dissipation and fluctuations even at zero temperature. 

Using a Runge-Kutta scheme, Eqs.[5] and [5] are solved 
sclf-consistently. In line with the discussion above, we 
find that switching occurs only for a limited range of 
mutually dependent parameters. The results in fig. [1] fo- 
cus on the power dependence of the magnetic response 
for uji = 10 15 Hz, uj = O.Slui, T = 0.2, and A = 10 4 Oe. 
While the response to laser fields Bi < 10~ 5 Hb is negli- 
gible on a femtosecond time scale (cl), it becomes signif- 
icant for 5 • 10~ 4 Hb (c2) and ultrafast reversal is found 
at higher fields of 10~ 3 H B and 2 • 10~ 3 H B (c3 and c4, re- 
spectively). At the highest fields, M z shows an almost in- 
stantaneous reversal accompanied by strong oscillations. 
The period of these oscillations is governed mostly by Hb 
while their decay depends on T. The reversal slows down 
with time and the approach to equilibrium depends on 
the coupling X(E) of the spins to the mode q which is 
proportional to the power. The curves c3- and c3n show 
the strong dependence of the reversal on the chirality of 
the laser. 
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FIG. 2: Switching paths for our theory (solid curve), the cor- 
responding modified Bloch equation MB, MB2 for l/|wj|T = 
0.06, and the Landau-Lifshitz (LL) equation for a = 0.25. 
Zero anisotropy has been assumed for MB, MB2 and LL and 
all other parameters are the same as in Fig. []Jc3. With- 
out anisotropy, the MB equations lead to the trivial solution 
M = 0. Dynamics of the transverse x-component of the mag- 
netization as predicted by our theory (dashed line), MB2, and 
LL are shown in the inset. 



Next we compare our results to those of the modified 
Bloch equations (MB, MB 2) 



dM 
~df 



1 



7 Mx (H + H B )--(M- X (H + H B )), (13) 
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and the (modified) Landau-Lifshitz equations (LL) 

— = 7 M x (H + Hg) - c* 7 M x (M x (H + H B )) . 

at 

(14) 

All equations are written in the rotating frame and will 
be solved without taking account of the anisotropy. The 
anisotropy is important only after the reversal of the 
magnetization and the laser is off. The common modified 
Bloch equations (MB) do not include the Barnett field in 
the relaxation term and have been successfully used to 
interpret low frequency free-induction decay experiments 
in paramagnets. However, as shown in Fig. [51 MB does 
not predict switching when the the Barnett field is much 
larger than the laser intensity. In contrast, the modified 
Bloch equations with the Barnett field accounted for in 
the relaxation, MB2, agree better with our results. 

The LL equations with the Barnett field in the relax- 
ation also give rise to fast reversal (Fig. [21 discontinuous 
line). The choice of parameters are made such that all 
solutions cross the point S at the same time. In the lab 
frame, the LL equations, unlike MB, also give reversal 
but for much longer times. For a — 0.25, a switching 
time of about 40 ps is found in the lab frame (not shown) . 
It can be shown, with little effort, that the Larmor equa- 
tion in the lab frame has a solution with a non-zero static 
z-component such that M z = cos 6 (cf. Eq. [TOf . 

The inset in figure [2] shows the relaxation of the trans- 
verse components of the magnetization. At short times 
both MB2 and LL equations fail to capture the real dy- 
namics of the reversal process. The failure of LL equa- 
tions to adequately describe the dynamics in our system 
clearly invalidates the assumption of the universality of 
the LL damping term as was claimed by Koopmans et al 
[l9j |. Note, that the LL relaxation term is also invalid at 
short times when the relaxation is caused by momentum 
relaxation of the conduction electrons [13] • 

Laser powers of the order of 10 15 W/cm 2 are relevant 
only to plasmas since most solids are ionized by such in- 
tense radiation [2(j. In this case, we have to abandon 
the identification of the mode q as that of an optical 
phonon mode and simply assume that it is a term due 
to collisions between the plasma particles. At these in- 
tensities, the spin exchange term S z a z starts to become 
comparable to the Barnett field. For powers in the range 
10 17 — 10 20 W '/ cm 2 , full reversal of the polarization of the 
magnetized plasma is not possible since the Zeeman en- 
ergy term grows much faster than the exchange term. At 
powers higher than 10 20 W/cm 2 , the exchange term dom- 
inates all other terms and full polarization of the plasma 
along the z-axis becomes again possible. 

Figure [3] shows that in the laser configuration studied 
here (k || z), the Barnett effect behaves similar to the spin 
momentum exchange term. The optical Barnett effect is 
therefore a spin momentum transfer effect and is not or- 
bital in character as in the original experiment of Barnett 
since it does not require charged particles. This can also 
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FIG. 3: Large field solutions for lui = — 10 15 Hz and Bi = 
5.0 Ht leading to an exchange field of 1.73-ffg . Initially the M y 
(middle) and the M z (lower)-components are zero. Solutions 
of the MB2 equations (dotted lines), with l/|u>j|T = 0.06, are 
plotted on top of our solution. 

be easily seen if we write the energy, Hi = Bf /An, of the 
laser wave as a Zeeman term, Hi — I/7V • J. The 
energy of the laser is not modified in the system studied 
here. Hence, the laser plays the role of a second bath 
that is strongly coupled to the spins. The spin momen- 
tum is transferred from the laser to the spin through the 
Barnett field as can be seen from the last term in Eq. [TJ] 
In summary, we have shown that circularly polarized 
light can induce femtosecond magnetization reversal, via 
spin momentum transfer, when optical phonon modes 
with frequencies comparable to those of the light are 
present. This fast reversal can be recovered qualitatively 
from MB and LL equations only if their relaxation terms 
are modified to account for the Barnett field. The switch- 
ing is found to be insensitive to the anisotropy which 
makes the proposed mechanism very attractive to high 
density magnetic recording. In high energy plasmas, an 
additional spin exchange transfer term between the laser 
and the polarization becomes dominant and can be used 
to control the induced magnetization. 
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